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Here I discuss a modified gravity tiieory. Modified gravity is often mentioned today in re- 
lation witii dark energy in cosmology. But I concentrate on the first systematic modified gravity 
approach to the mass discrepancy in galaxies and clusters, Milgrom's modified Newtonian dy- 
namics (MOND) [|, I, |l, and report on a relativistic implementation of it [Q]. My approach is 
phenomenologically motivated, not theory motivated. 



1. Introduction 



There is missing mass in disk galaxies. We see this from the rotation curve, the plot of the 
circular orbital velocity with radius. A few kiloparsecs out from center of a disk galaxy, the rotation 
curve flattens out. From the observed stellar distribution Newtonian theory would rather predict a 
falloff. Interpreting the rotation curve a la Newton tells us that away from the central regions, the 
mass enclosed in a given radius grows proportionately to radius, ultimately giving a galaxy mass 
an order of magnitude higher than the visible one. It is as if some invisible matter contributes to 
the galaxy's gravitation; hence the dark matter (DM) hypothesis. This finding is accompanied by 
a mysterious correlation between infrared luminosity of a disk galaxy Lk and the asymptotic rota- 
tional velocity v^: Lk°^ Vq^ (Tully-Fisher law) [||, |6|]. The scatter of the data about the Tully-Fisher 
law is consistent with observational uncertainties, i.e., there is no intrinsic scatter — a wonderful 
situation in astrophysics, and one hard to understand from the point of view of the DM paradigm. 

And several lines of evidence show there is missing mass in clusters of galaxies. Random 
galaxy motions in a cluster define a (dynamical) velocity scale. The virial theorem then gives an 
estimate for the mass interior to the cluster, which comes out much larger than the sum of galaxy 
masses [Q, §]. In the last two decades hot X-ray emitting gas has been found in abundance in 
many clusters. From its temperature distribution not only the cluster mass but its distribution can 
be inferred. The extra mass in gas has only ameliorated, not removed, the missing mass problem: 
there is still a factor of ~ 5 discrepancy [^]. Some clusters are seen to lense distant background 
galaxies gravitationally. The measured light deflection has been used to determine the mass and 
its distribution in each cluster. Again, the mass comes out larger than the visible mass. When the 
dynamical and lensing methods can be compared, they give a similar mass discrepancy factor. In 
certain circles this is considered a vindication of the DM scenario. 

The DM paradigm [||], now 30 years old, holds that each galaxy is nested inside a dark halo 
whose slowly declining mass density profile makes the rotation curve approximately flat over some 
range. An halo model necessarily includes at least three free parameters. They may be the ve- 
locity dispersion of the dark particles, the inner cutoff radius (indispensable to prevent mass di- 
vergence), and the mass to luminosity ratio M/L of the visible matter. But dark clouds hover 
over this paradigm. For example, fine tuning of the halo and visible galaxy model parameters is 
needed to avoid a never-observed cusp in the rotation curve ||^. And attempts to understand the 
Tully-Fisher law as reflecting galaxy genesis [ IC ] have problems in explaining how such inherently 
messy processes can give a sharp correlation [11]. Add to this the embarrassment of the giant 
elliptical galaxies; now that extended rotation curves have been determined for some of these, the 
data reveal no evidence for a dark halo |p^]. And yet it is widely accepted that elliptical galaxies 
are formed by merger of disk galaxies. One wonders, where the DM is, here or elsewhere ? No 
experimental or observational claim of direct discovery of "dark particles" has been sustained. 
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2. The MOND paradigm 

The MOND paradigm is two decades old [[I], ^ ^. Before MOND many workers toyed with 
the idea that a change in the l/r^ law at large length scales is the explanation for the mass discrep- 



ancies []13[]. But as Milgrom first realized, any modification attached to a length scale would cause 
the larger galaxies to exhibit the larger discrepancy [|l|, [Tl]]. This is contrary to the observations: 
there are small low surface brightness galaxies (LSBs) with large discrepancies, and very large 



spiral galaxies with small discrepancies [14]. 

Milgrom suggested acceleration scale as the more relevant one. Since then the presence of 
an acceleration scale in the observations of the discrepancy in spiral galaxies has become more 



evident [11, 14|. For example, consider the log-log plot of the dynamical M/Lk vs. the radius R 
at the last measured point of the rotation curve for a uniform sample of spiral galaxies in the Ursa 
Major cluster [^. (Assuming a spherical mass distribution, one calculates M with the Newtonian 
formula M = Va^R/G). The result is a scatter plot. But if the abscissa is changed to centripetal 
acceleration (a = Va^ /R ) at the last measured point, there is a nice correlation: M/Lk °^ 1/a for 
a < lO^^cms^^ and M/Ljc « 1 for a > 10^*^cms^^. Now population synthesis models suggest 
that M/Lk should be about unity [|l^]. So it is clear that the Newtonian analysis fails for a > 
lO^^cms"^. 

MOND maintains that the familiar equality between a test particles's acceleration a and the 
ambient Newtonian gravitational field — VO/v is not exact, and is replaced by 

^(|a|/ao)a = -V<I>;v (2.1) 

where oq is an acceleration scale of order 10^*^cms^^, and /^(x) x for x ^ 1 and asymptotes to 
unity for X » 1. So in laboratory settings, where \a\ ^ ao, = I and we get back to Newtonian 
physics. 

Now, when |a| <C ao, which occurs frequently in galaxy outskirts, we have the extreme MOND 
equation lala/ao = — V<I>Ar. For circular orbits on a plane we can set |a| = jr; both sides of the 
equation then behave as 1/r^ provided v = const, (flat rotation curve). Equality of the coefficients 
gives M = (Gao)^v^, or Lk = {Ga^M /Lk)^^v'^ . On the simple assumption, which accords well 
with the population synthesis models, that M/Lk is constant, this is precisely the TuUy -Fisher law. 
Notice that the TuUy-Fisher law is derived as a sharp relation, just as observed. 

Not only this but MOND has been successful in fitting, for over 10^ spiral galaxies of all 
sizes, the detailed shape of the rotation curves using the observed distribution of stars and gas 
as inputs [^; in these fits the only free parameter is M/Lk (or if more convenient M/Lg in the 
blue band). For 90% of the said galaxies there are no significant glitches, and the best fit values for 
M/Lk or M/ Lb are in excellent agreement with those predicted by the population synthesis models. 
Comparable success in fitting rotation curves can be had for dark halo models only by adjusting 
the above mentioned three halo model free parameters. MOND is definitely a more parsimonious 
hypothesis than dark halos. 

We mentioned that the dynamics of some giant elliptical galaxies give no evidence for dark 
matter. In MOND this makes sense: these galaxies are found to have accelerations that exceed ao 
well out of the core. To first approximation there should thus not be a deviation of MOND from 
Newtonian dynamics in them. Milgrom and Sanders have taken a more precise look at the giant 
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elliptical NGC 3379 (dynamical data from 2003) and shown that it is well described by a MOND 
based model of 1983 vintage [16|. 

But it is not all roses; there are some embarrassments for MOND. Let us mention three here [|l|, 



17| , |18| ]. A system of point masses moving a la MOND does not conserve momentum or angular 



momentum. If that were not enough, the motion of a star orbiting circularly in the outskirts of 
the the galaxy is predicted to be (extreme) MONDian if the star is regarded as a point mass, but 
Newtonian if attention is paid to the high accelerations to which the star's gas parcels are subject 
to. Finally, open star clusters in the Galaxy, e.g., the Pleiades, are all characterized by internal 
accelerations small on scale ao, yet display no mass discrepancy whatever, quite at variance with 
MOND's intent. 



3. Nonrelativistic field theory for MOND 

In light of the above, what are we to make of MOND ? A sector of the astrophysical community 
regards it as just a good phenomenological summary of galaxy phenomenology. But it is clearly 
more interesting and promising to consider MOND as a reflection of unorthodox physics. Here 
we shall consider how far it can be due to modifications of standard gravity theory. (Alternatively, 
Milgrom has considered MOND as a modification of Newton's second law ^ [T^, but it is unclear 
to me how such scheme would translate to the relativistic domain). 

The first step, a modified nonrelativistic gravity theory, dubbed AQUAdratic Lagrangian the- 
ory (AQUAL), solves the mentioned three problems with MOND. AQUAL is based on the la- 



grangian [18] 



/(Ll)-p* (3^1) 



where is the true gravitational potential (meaning that for a test particle a = — VO), and p denotes 
the total matter mass density. The very fact that the formulation is lagrangian based automatically 
removes the problems with conservation laws. 

The root of those problems becomes clear when one derives the equation for <I>: 

V-[a(|V<I>|/ao)VO]=47lGp; (3.2) 
HVy)^dfiy)/dy. (3.3) 

Comparing this with the Poisson equation 

V • VOa, = 47iGp (3.4) 
we see that //(|V<I>|/fl'o)V<I> = V^n — V x h, with h a vector field which insures that both sides of 



the Eq. Eq. ( |3.2[ ) have the same (nonvanishing) curl. Replacing — > —a we have the AQUAL 
equation 

//(|a|/ao)a = -V<I>N + Vxh. (3.5) 

With the proviso that /)(x) and ao here are the same as in Milgrom's scheme, this is a corrected 
version of the MOND equation (p]). The correction V x h is what brings the dynamics into 



harmony with the conservation laws, as verified by direct calculation [18]. 
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The same type of calculations show that the center of mass of a collection of mass points 
moves in the averaged ambient — V<I> in accordance with Eq. Eq. (^^), regardless of how strong the 
internal fields acting on the individual mass points are. This brings agreement with the observation 
that dense stars and tenuous gas clouds seem to trace out exactly the same rotation curve for each 
galaxy studied using both. The AQUAL equation, unlike the original MOND equation Eq. (^), 
respects the weak equivalence principle (all objects move the same way on a given gravitational 
field). 

The open clusters paradox is also solved by AQUAL in harmony with Milgrom's early con- 
jecture [Q] that the environment of a system can suppress its MOND-like behavior. Calculation 
reveals that the dynamics of a subsystem with internal accelerations below oq will approach New- 
tonian dynamics if the enveloping system subjects it to a — of magnitude exceeding ao [p^]. 
This is precisely the situation for the open clusters, many of which are embedded in the gravita- 
tional field of the Galaxy where it is not weak on scale ag. They should thus behave in Newtonian 
fashion, as found. 

Further, it is found that if both the internal accelerations and the stronger external field gg = 
I V<I>| acting on a subsystem are weak on scale ao, the subsystem's dynamics is again Newtonian, but 
with an effective gravitational constant G/fi{ge/ao) [|l8|]. Now the dwarf spheroidal companions to 
the Galaxy, Leo II, Sextans, Draco, etc. have low internal accelerations, and are so far out in the 
Galaxy's outskirts that they feel a weak external field. The prediction is then that they should have 
Newtonian dynamics (which they seem to do in the sense that the dynamically inferred mass dis- 
tribution mirrors the luminosity distribution), but the M/L determined by virial theorem estimation 
should be larger than expected from the stellar population involved by a factor \/^{ge/aQ). And 
indeed, the dwarf spheroidal exhibit anomalously high virial M /L\, a fact anticipated by Milgrom 
in his classic papers [^. 

For spherical symmetry the vector h vanishes; hence for spherical galaxies or clusters the 



AQUAL formula Eq. ( [3.5| ) reduces exactly to the MOND formula Eq. ( |2. 1| ). We have already 
mentioned that MOND agrees with the recent finding that the dynamics in giant elliptical galaxies, 
many of which are quasispherical, do not require dark matter for their explanation when interpreted 
by Newtonian theory. We may also recall an early MOND prediction [^] that spherical stellar 
systems with equilibrium isothermal velocity distributions can exist only with typical accelerations 
« ao or below. To be sure this is found to be true for the giant molecular clouds in the Galaxy, its 
dwarf spheroidal satellites and a fair majority of X-ray bright clusters of galaxies. Giant ellipticals 
and globular clusters can have higher accelerations, but it is known that the velocity distributions 
here are anisotropic. It is a striking finding that all the five classes of objects have central 
accelerations within a factor of five of Milgrom's constant ao. originally defined with respect to 
disk galaxies. 

And what about disk galaxies, MOND's original playground ? It was known early that the h 
term in Eq. (3^) falls off faster with increasing radius than the term. Numerical calculations 
by Milgrom show that in the inner regions h contributes only some 10 — 15% of the accelera- 
tion [21 1. This means that the old-fashioned MOND reconstruction of rotation curves from the 
visible matter distribution, which is known to be successful, is very much what a full AQUAL 
treatment would give. 

To be sure MOND and AQUAL have a common failing. Great clusters of galaxies, which 
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display accelerations up to a few times ao, show an average factor of five discrepancy between 
the Newtonian dynamical mass and the mass in stars plus gas. With MOND the discrepancy can 
go down to a factor of two, but does not go away. This has to do with the large accelerations 
in clusters |^]. AQUAL does not solve this problem because many of the clusters in question 
are spherical, making MOND fully equivalent to AQUAL. And just hke MOND, AQUAL does 
not give a prescription for calculating gravitational lensing by extragalactic objects, so that the 
mass discrepancies evidenced by the lensing of distant galaxies by clusters cannot be examined 
critically. Nevertheless it is clear that in most cases AQUAL theory does well phenomenologically, 
while having conceptual advantages over the pure MOND prescription. 



4. Steps towards a relativistic modified gravity 

Once MOND's philosophy is accepted, a relativistic implementation of it is indispensable, 
not only for the computation of gravitational lensing, but also to permit investigation of other 
relativistic systems like close binary pulsars or cosmology. I now show how hints from AQUAL 
led to the construction of a good candidate for a relativistic gravitational field theory for MOND [Q]. 

Such theory cannot arise from replacing the Einstein-Hilbert Lagrangian density R of gen- 
eral relativity (GR) by some f{R). True, part of R is quadratic in derivatives of the gravitational 
potentials (metric components), and such replacement is suggestive of a nonrelativistic limit like 
Eq. ( |3.l[ ). But R also contains a term involving second derivatives of metric components. In the 
Einstein-Hilbert action this term reduces to a boundary term and is irrelevant. But in f{R) it would 
not do so, and would actually dominate the previously mentioned term. To the above "poor man's 



argument" one can add the Soussa-Woodard theorem [ |22| ] that a gravitational theory built on metric 
alone, whether with global or local lagrangian density, cannot both have a MOND like nonrelativis- 
tic limit and reproduce the anomalously large gravitational lensing seen in clusters of galaxies. 
These problems lend credence to the first relativistic approach to MOND, relativistic AQUAL 



(RAQUAL) theory [18, 23]. It promotes <t> in AQUAL to a scalar field (|), and makes its lagrangian 



density (henceforth I set c = 1 , but display G) 




where is the inverse of the Einstein metric g^p and / may be different from /. One adds to 
this the Einstein-Hilbert lagrangian density R as well as a matter lagrangian density L,„. This last is 
obtained from the usual curved spacetime one by the replacement g^p —i- g^p ^ e^'^gaf,- addition, 
when integrating £,„ over spacetime to form the matter action S,n, we do it with the determinant of 
g(xp. We call goip the physical metric. For example, a point mass m has the action 

S,n = -m J (-g„p J^«J^P)i/2^ (4.2) 

obtained by replacing g^p — > g^p = e^^g^p in the familiar integral over the worldline. The field 
equation deriving from / Li^{—gy/^d'^x + Sm is [p{y) = f{y)] 

/ _aBjv J. \ 1 

47iGg«Pf„p (4.3) 
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where T^p is the energy-momentum tensor obtained by varying the matter action with respect to 
gg.p. For correspondence with AQUAL we choose f{y) = for < j ^ 1 and f{y) y for 
j> 1. 

How does RAQUAL yield Newtonian and MOND limits where relevant ? As well known, in 
the nonrelativistic limit slow motion in GR is governed by the temporal-temporal metric component 
which can be approximated by — (1 + l^ii), with coming from the usual Poisson equation 
with p as source. According to Eq. g,;'s job is taken over by f,, — (1 + + 2(|)) so that it 
is O = + (|) which represents the nonrelativistic gravitational potential. Thus (|) will replace the 
gravitational potential of the dark matter. For quasistatic situations with nearly flat metric, its field 



equation has the form of AQUAL's Eq. (|3.2D. 

Now suppose that because of the large scale of the system, |V(|)| ^ qq. Comparison of (j)'s 
equation with Poisson's Eq. (3^) shows, in light of the fact that here /J <C 1, that |V(|)| S> |V<I>a?|. 
Consequently V<I> « V(|), so that the physical potential approximately obeys the AQUAL equation 
Eq. This means that for weak fields (MOND limit), RAQUAL reduces to AQUAL, and 

inherits the latter 's successes. Now when the system is such that |V(|)| S> we shall have « 1 
meaning that ([) approximately obeys the equation for <I>a?. Thus V<I> 2V<I>/v. This just mean that 
at the nonrelativistic level we have Newtonian gravity, but with an effective gravitational constant 
Geff = 2G. We shall not dwell on whether this result is also consistent at the post-Newtonian level 
because a tenacious problem renders RAQUAL unviable. 

Because f „p = e^^g^^}^ (conformally related metrics), all optical phenomena, light ray paths 
included, are exactly alike in the two metrics. Any hope for anomalously strong light bending, as 
required by the gravitationally lensed galaxies, hangs on the Einstein metric being substantially 
different from that in GR with no dark matter. However, in the case of the clusters |V(t)| = o(ao)> 
so that the scalar's energy density is of order G^^a^^ ~ G^^Hq^ (here Hq is the present Hubble 
"constant"). This cosmological-scale energy density cannot compete with that furnished by the 
stellar component in a cluster. Hence, in RAQUAL ^oip is not very different from that in GR, and 
no anomalously strong lensing can be expected from it. In brief, RAQUAL fails the lensing test 
miserably. 

Sanders and I tried to rescue the RAQUAL idea by replacing the conformal relation between 
metrics by g„p =exp(2(|)) [gap + fn(|),a(|),p], where (D is some function of the invariant g"P(|)_o[<t>.p [@, 
25]. Here gap will bend light rays differently from g^p. We only assumed that U5 <0, and this 
to guarantee that g^^p and g^cp both have the usual signature, and that weak gravitational waves 
(perturbations of gap) propagate inside or on the physical light cone (that of gap). A calculation 
then showed that if a cluster can be approximated as static with d(^/dt = 0, then the light bending 
is weaker than in GR for the same mass: the lensing problem becomes worse than in RAQUAL ! 

An alternative proposal by Sanders [26] has been more fruitful: in AQUAL multiply just the 
part of the metric orthogonal to a fixed unit time-directed vector il", namely g^p +itailp, by e^'^'^ 
and that part coUineai" with il", namely ilaiip, by e^'^, taking the sum of the two parts as the physical 
metric. In this so called stratified theory, as in a little known pre-GR paper by Einstein, the ruse 
immediately produces extra light bending. Indeed it can resolve the anomalously large lensing 
problem. But ita defines a preferred frame, thus breaking local Lorentz invariance, and because 
ila is specified as time-directed in all situations, the theory cannot be generally covariant. The 
first problem is not critical as ila couples to matter only through the physical metric, so that the 
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breakdown of local Lorentz invariance is restricted to the gravitational sector. Sanders found that 
the breakdown is only weakly reflected in the relevant parametrized post-Newtonian formalism 
coefficients. The lack of general covariance is more serious. This problem cannot be resolved until 
a prescription for determining ila is given which makes no reference to particular coordinates. This 
is where TeYe S takes over. 

5. Structure of the tensor- vector-scalar modified gravity 

This covariant modified gravity theory, dubbed TeYeS for short, starts off from Sanders' pre- 
scription for the physical metric: 

fap = e-^^{gap+ilailp) " e^^Hiailp (5.1) 

where one demands that g^Pilailp = — 1 , but does not stipulate any particular orientation for the 
vector. The Einstein metric gets its dynamics from the usual Einstein-Hilbert action Sg. The matter 
action 5,,, is built by replacing the Minkowsky metric and various derivatives by the physical metric 
and covariant derivatives with respect to it. 

The vector field's dynamics come from the postulated action 

5v = / [g«Pg^^%,,]H[p,v] - 2(V^:)(g^^il,iiv + 1)] {-gY^^d'x; (5.2) 

the normalization of ii^ is enforced via the Lagrange constraint with multiplier function X. The 
is a dimensionless positive coupling constant of the theory. Note that it is the Einstein metric 
which appears in 5y. To derive the vector's equation one is directed to vary it^, as opposed to the 
contravariant components (11^ = g'^^Uy), in 5,, + 5',,,. Now ily enters 5,,, through the where 

= eV + 2rtP%ilvsinh(2(^). (5.3) 

This and {—gY^^ = e^'^'^{—gy^^ have to be varied giving thereby source terms depending on the 
energy-momentum tensor [(— g)'^^7kp = — 28Sm/8g"P]. The resulting equation, 

isTill^'P] ;p + XU« + 87cGa2ilP(^,p^«^(^,y = 87iG(l - e-^^)g''^'ii%p, (5.4) 

is best regarded as giving the value of X together with those of three components of il^ (the fourth 
coming from the normalization condition). 

Interesting consequences of Eq. ( ^^ surface if we model matter by a perfect fluid: 

Tap = P"aMp +P{gaf, + "aMp) (5-5) 

Here p is the proper energy density, p the pressure and Ua the 4-velocity, all three expressed in the 
physical metric (m" = g^^wp). It is found that both for a static system, e.g. a quiescent galaxy, and 
for isotropic cosmology, is collinear with the fluid velocity m"; in particular, they both point in 
the time direction. To put it picturesquely, If^ "seeks out" the frame in which matter is at rest. In 
this manner is Sanders' desideratum obtained dynamically rather than by fiat. 
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The scalar's action is taken of the form 



(5.6) 
(5.7) 



where k is another dimensionless positive parameter, ^ is a scale of length, and F is some function. 
The tensor is used instead of by itself in order to prevent causality problems (see Sec. 
The a in the action is a nondynamical auxiliary field which allows the kinetic term of (|) to be 
presented as quadratic, which it really is not. This a is determined by the extremum condition 
55v/5a = 0. The (j) equation is obtained by varying Ss + Sm with respect to (j). Now (|) enters Sm 
through the g^^, and variation of this last introduces source terms containing the energy-momentum 
tensor. In terms of the function ^{y) (to be distinguished from Milgrom's p) defined by 



-luF{ju)--f?F'ifx)=y, 



we obtain the equation 



(5.8) 



(5.9) 



which is plainly reminiscent of the AQUAL equation Eq. (3.2). 

The equations for the metric result from varying in Sg + 5„, + 5v + 5^. We observe that Sm 
is built from part of which varies in accordance with Eq. (5.3 ), and thus introduces the matter's 
energy-momentum tensor into the equations. These are 



G, 



SnG 



Tap + (1 - e-4^)H^f^(„%) + t„p1 + 0„p 



(5.10) 



where 



?ap- 



(5.11) 



It is most convenient to vary rather than g^p. If the latter procedure is implemented cav- 
aUerly, it is easy to get equations which are not the raised indices versions of those coming from 



the former procedure. C. Skordis has observed that this happens because in writing Eq. ( |5.3[ ) as the 
inverse of Eq. (5J), one is already assuming that g^Pilailjj = — 1 before performing the variation, 
whereas this normalization should be left free to be enforced by the Lagrange multiplier term. At 
any rate, the difference in the equations obtained boils down to a difference in the X appearing in 



the vector's 0(xp» but since X has to be computed after the fact from the vector's equation Eq. (pl4|), 
and this last is also changed, the difference disappears in the final equations. 

It will be observed that F shows up in the scalar's equation and contributes to the source of 
Einstein's equations. In my original paper I made the choice 

3 



-2 



which leads to 



H {A + 2n-A^? + ^?)+2\n[{l-nY] 



-2)2(1 -^)- 



(5.12) 



(5.13) 
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The important features in this y{iuj are that for small y, /uoc and that large y occurs for ^ 1. 
There are other choices of F (/li) that will give these, so the above choice is to be regarded as a toy 
model. 

Other opportunities for generalizing TeVe S exist. For example, one can use f in place of g""^ 
to construct 5,, and 5^ . Or one can replace the coefficient of in by any number below — 1. 
We have not yet elucidated how the theory's predictions would be changed then. In addition, one 
can contemplate kinetic terms in Sy containing the symmetric combination ii[a;/.i)^{^;v) opposed 
to the antisymmetric ones. But if such a term is added, as in the Jacobson-Mattingly modified 
gravity with preferred frames [27], the gravitational equations then contain second derivatives of iia 
coming from covariant derivatives; these complicate the formulation of the initial value problem. 



Finally one could add a kinetic term for a, as in the phase coupled gravity (PCG) theory [23 



28]. However, PCG tends to induce instabilities in bound systems [p9|], so this addition may be 



inadvisable. 



6. Nonrelativistic limit of TeYe S 

One reason for believing that TeYe S is not a pathological theory is that it has GR as a limit. 
Milgrom has remarked that by suitably rescaling (|) and a, one can see that as i oo and K 0, 
the action Sg + Sm + + reduces to its first two terms, with X vanishing and the two metrics 
becoming identical. This is pure GR. Physically this means that for systems small as compared to 
£ and for small K, TeYe S should be a passable approximation to standard gravity theory. We thus 
suspect that in the real world K <^ I and the scale i is cosmologically large. A similar comment 
applies to a second way to get the GR limit (which is established as yet only for static systems and 
isotropic cosmology): take k —> with i oc ^^^/^ and K oc k. 



In the nonrelativistic approximation of TeYeS, one first uses the vector equation Eq. (5.4) to 
compute X; one then substitutes it in the linearized Einstein equations, and discards all of Tap and 
©oip (apart from the X term) because they are quadratic in presumably small quantities like (|) « or 
p (the F term in T„p also turns out to be negligible after more detailed analysis [Q]). In terms of 

E = e-^^'{l+K/2)-\ (6.1) 

where <^c is the asymptotic boundary value of (|) (set by the cosmological model in which the system 
in question is embedded), one gets by the familiar procedure 

gn--il+2E^N)- (6.2) 

The coefficient of is not 2 as in GR because T^p appears in the source of Einstein's equations 
in three distinct forms, only one of which is the familiar one from GR. After multiplication of gtt 
by e^^ ^1+2^, 

f„ w -(1+20); <I> = S<I>;v + ^. (6.3) 

The physical gravitational potential here is quite close to that in AQUAL; recall that we expect 
K to be small, and it can be shown that (|)c is small in suitable cosmological models. Hence S is 
close to unity. 
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Let us now compute (j) for a static spherically symmetric system with the matter represented 
by the fluid energy-momentum tensor Eq. (^). Taking gg,p as flat, the scalar equation Eq. (| 
reduces to 



^^(ke^{V(^f^V(^ =kGp. 



Comparing with Poisson's equation Eq. (3.4) we see that 

V(|) = (/c/47i^)V<I>A.. 



We may now calculate from Eq. ( p3[ ) and cast the result as 



(6.4) 



(6.5) 



(6.6) 



The first equation is precisely the MOND equation (2.1) since by the nonrelativistic limit of the 



geodesic equation, and the form of gtt in Eq. ( |6.3| ), a = — VO. Actually in light of Eq. ( |6.4| ) /j is 
regarded as a function of |V(|)|. However, using the last two equations we may reexpress it, as well 
as /i, as functions of | V<I>|, just as required by MOND. 

Let us now examine the case where because of the weakness of the gravitational field, ^ ^ 1 . 
According to Eq. ( |5.13 ) the auxiliary variable y, the argument of fi in Eq. (5^), is y « 3fp-. Thus 
fi ^ {k/3y/^£\V^\. Using this in Eq. ( |63| ) and Eq. ( |6^ to eliminate VO^r and V(|) as well as ju gives 



A/=(A:/87lE)(-l + Vl+4|V<I>|/, 



ao = (3/S:)^/2(47lS 



(6.7) 



The positive root was chosen to get /j > 0. Of course, this expression for /j can be believed only 
when it gives ^ 1; assuming k is crudely o(l), this will be true when |V<I>| <^ao,m which case 
/J (^/47iS)|V<I>|/ao. But then Eq. (|6^ informs us that jj 'E^^\V<t>\/aQ. In fact, we get 



(6.8) 



|V<I>|V<I>/ao = SV<I>^ 



If one identifies in Eq. ( p77| ) with Milgrom's constant, this is precisely the extreme (|a| ^ ao) 
limit of the MOND equations (^?1|). The apparently superfluous factor S is understood by referring 



to Eq. ( p.2\ ) which shows that in TeYe S it is SOat which plays the role of nonrelativistic potential 
in the Einstein metric. More on this below. 

Consider now the case ^ « 1 which means y = V(t)|2 > 1. It follows from Eq. that ji w 
{E + k/4n)^^. And that same equation tells us that = {E-\-k/4n)V^N. The proportionality of 
the two gradients means that the acceleration —V<t> has Newtonian form, so that ^ 1 corresponds 
to the Newtonian limit of TeYeS. And from the coefficient of V^n we learn that the physical 
gravitational constant in the Newtonian regime is Gn = (S + k/4n)G. The contribution EG has 
been mentioned already as coming from the gravitational equations; the part kG/4n is contributed 
by the scalar field. 

How big are the departures from exact Newtonian gravity in TeYe S near = 1 ? Let us expand 
y(ju) about the singular point = 1 in a Laurent series: 



3;=_(l_^)-l + 0(l_^). 

4 



(6.9) 
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Because there is no constant term here, ^ = 1 — ^y^^ = 1 — V(|)p)^^ furnishes an accurate 
inversion. Eliminating in favor of |V<I>| with help of Eq. ( |6!^ (where we set ^ = 1) and 
Eq. gives 

Ai^T^ 1-^^ • (6.10) 



In interpreting this result one must realize that if in all formulae involving observables, the 
constant G is reexpressed in terms of G^, the prefactor G/G^ here will disappear. Thus the leading 
term in Ji is really unity, as in MOND. Next we have a correction which still relates to nonrelativistic 
gravity, and vanishes as |V<I>|/fl;o ^ °°. In a strong gravity system, i.e., the solar system, this 
correction is small enough not to bring about conflict with the very accurate ephemeresis ii k = 
O (1). But a much smaller k may cause trouble. So far the only reason for being interested in small 
k is cosmology, where small k leads to easy to understand models. But if good models cannot be 
had for ^ = o(l), one would have to rethink the form of F{ij) and perhaps other features of the 
theory. 

The cases ^ and 1 conform with the pure MOND equation (2.1). But when ^ has 
an intermediate value, TeWe S introduces corrections to MOND. In this regime <I>a? and (j) make 
comparable contributions to <I>. But they satisfy two different equations, (3.2) and (3.4), which 
means that in general is not exactly a solution of a MOND type equation. 

As I mentioned in Sec. ^, there are clusters of galaxies whose dynamics are not well described 
by MOND. One possible reason may be that such clusters harbor still unseen matter [pO|]. But the 
above mentioned correction offers an alternative. Even for spherical clusters TeVe S predicts devi- 
ations from pure MOND behavior, particularly in those systems with typical internal accelerations 
near a^. This possibility has not yet been explored quantitatively. 



7. Relativistic corrections in Te\e S 

The post-Newtonian effects have often been used to rule out competitors of GR; how does 
TeNe S fare in this respect ? Here I briefly sketch how to calculate from TeVe S the parametrized 
post-Newtonian coefficients (3 and y for the physical metric exterior to the Sun, 

ga^d/'dx^ = -eVf^ + e^[Jp2 + p2(je^ + sin2eJcp2)], (7.1) 
e"^ = \ - IGNinp-^ +l^GN^n?p-^ + 0{p-^), (7.2) 
= \ + 2yGNmp-^ + 0{p-^). (7.3) 

Here m is the Sun's physical mass and G/v was defined previously. The y determines the light 
deflection and radar signals time delay, while both it and (3 determine the perihelion precessions of 
the various planets. 



The first stage is to write the Einstein metric in analogy with Eq. ( [7.1D with the metric functions 
expanded as 

= l-r,/p + a2(r,/p)2 + .... (7.4) 
e^=\ + ^,rjp + ^2{rg/9f + --- (7.5) 



012/12 



Modified Gravity vs Dark Matter: Relativistc theory for MOND 



Jacob D. Bekenstein 



The Tg is a scale of length to be determined, and so one coefficient can be taken as — 1 (negative 
because gravity is attractive). Symmetry and normalization dictate that 

H« = {e^/2, 0,0,0}. (7.6) 

With this one solves the scalar equation to get 

c^(p)=c^.-^ + 0(p-^). (7.7) 

Again, <^c is the asymptotic value of (|) and nig is a particular integral over p and p of the sun's matter 
which is close to the everyday formula for mass iQ]. 

Next one computes outside the Sun the parts of Tap and ©oip which fall off as p^^. Using them 



one solves the G/ and GpP components of the gravitational equations Eq. (5.10) by power series to 
determine the coefficients (with errors in Ref. [H] corrected) 



Pi = l; cx. = i; h = - + -K-'^. (7.8) 

^8 



2 8 16 87lr„2 



The next stage is to match the interior and exterior solutions for at the Sun's surface, which 
procedure determines that = IGm^ + O {rg'^ /R) where R is the Sun's radius and nig is an integral 
over the matter variables, T„ and 0,,, which agrees with the naive expression for mass in the Sun 
(as does m,) to 10^^ fractional precision. 

Finally one uses Eqs. {5A), (7^) and ( [7?7| ) to calculate the physical metric. The e^"^^^ factors 



entering in transformation Eq. (5.1) are gotten rid of by redefining units of length and time. The 



result can be put in form Eq. ( [7.1[ ) with the values (corrected for an error in Ref. [Q]) 

y=l; p = l, (7.9) 

where m = {G/GN){mg + kms/4n) is the observable mass. Both P and y are here the same as 
in GR. The classical solar systems tests cannot thus distinguish between the theories with current 
measurements precision, and TeYe S passes the elementary post-Newtonian tests. A desirable thing 
for the future is a calculation for TeYe S of the three post-Newtonian coefficients which have to do 
with preferred frame effects, and are today well constrained by experiment. 



8. Gravitational lensing according to TeVe S 

We recall that lensing by clusters of galaxies gives evidence that the same mass is responsible 
for galaxy dynamics in and gravitational lensing by the clusters. If there is DM with suitable 
distribution, this accords well with GR. As we now show, TeYeS does just as well [Q]. 

It is sufficient to work in linearized gravity. Write the metric as 

ga^=^ap + hap-'-^ap^'%b (8-1) 
where the perturbation from flat spacetime, hy^, satisfies the Lorentz gauge condition 

riP%/j^8 = 0. (8.2) 
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Then the Einstein tensor to O {h) takes the form 



Gap = --11^^9y9s/Zap- 



Accordingly, the gravitational equations (5.1C) take the form 



11^^ /%p = - 1671G [f„p + ( 1 - e-4^)ll^f^(„%) + x„p 



20, 



'ap 



(8.3) 



(8.4) 



Galaxies and clusters of galaxies are quiescent systems, so we may drop time derivatives of 
/lap; likewise, we drop the quadratic tensors Tap and ©ap apart from the term — AiXailp which is 



not ostensibly small. The X is calculated from Eq. (5.4) after neglecting the (|) terms. For Tap 



we take the form Eq. (5.5). Not surprisingly we recover an expression for hu we already met in 
Sec. htt = —AZ^^, where comes from Poisson's equation with p as source, and E has been 
defined earlier. Now in nonrelativistic systems ^ ^ p while velocities are small, so that the spatial- 



spatial and spatial-temporal components of the source of Eq. ( |8.4| ) are small compared to that of 
the temporal-temporal one. We can thus claim that hjj and htt ^ 0, so that 



5ap = (1 - 2S<I>;v)'nap -4S<I>;v5a5p. 



(8.5) 



We now use Eq. ( p.l| ) to linear order in ([) to construct the physical metric. The corresponding line 
element is 



where again = SOat + (j). 



-(1 + l^)df + (1 - 2<^)5ijdx'dxj 



(8.6) 



Metric Eq. (8.6) is the usual linearized metric used in GR to calculate extragalactic lensing. So 
just as in GR, in TeYeS a single potential is responsible both for galactic dynamics and for lensing. 
This solves the principal observational problem of early modified gravity theories in regard to 
lensing (Sec. ^. Of course within GR comes from the Poisson's equation with DM included in 
the source, whereas in TeVe S both the <I>Ar and (|) parts of <I> have only visible matter as sources, and 



the latter is computed from Eq. ( |3.2[ ) which is a suitable limit of Eq. ([5.9|). These two prescriptions 
for need not agree a priori, but as we have argued, nonrelativistic low acceleration dynamics in 
TeYeS are approximately of MOND form, and MOND's predictions have been found to agree with 
much of galaxy dynamics phenomenology. We thus expect TeVeS's predictions for gravitational 
lensing by galaxies (and some clusters of galaxies) to be in no way inferior to those of GR dark 



halo models. Now the MOND formula Eq. ([2. 1|), and TeYeS, at least with our choice of F, both 
predict that asymptotically the potential of an isolated galaxy grows logarithmically with distance 
indefinitely provided the environmental effect discussed in Sec. |3| is not important. Dark halo 
models do not. So TeYeS for a specific choice of F is in principle falsifiable. Dark matter is 
less falsifiable because of the essentially unlimited choice of halo models and choices of their free 
parameters. 



9. The problem of causality 

Not only did the RAQUAL theory of MOND fail on account of its lack of specific gravitational 
lensing, but it also exhibited superluminal scalar waves. In fact, modified gravity of any sort with 
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a scalar sector is easily afflicted by such acausality [29|. TeVeS, although designed for other 
purposes, resolves this problem neatly. 

To show the root of the problem and its solution, I will focus on the scalar equation, Eq. (^^), 
in a region outside the sources. Separating (|) background and perturbation as (|) = (|)b + 5(t), and 
linearizing in 5(|) gives 

g«p - + 2^//«//P) a„ap 5^ + • • • = 0, (9.1) 



where = {h^'^d ^fjd ^v^)^ ^^^h^-^d ^p^, ^ = dln^(y)/dlny, the • • • stand for terms with one (|) deriva- 
tive (for which reason we display only plain and not covariant derivatives), and the subscript B is 
dropped here and henceforth. If we work with RAQUAL's Eq. (p3|), the only difference is that the 
term -il«UP is missing in Eq. ( |9. 1[ ) and in the definition of In a static situation = 0, and 
as mentioned in Sec. ^, has only a time component; hence for both theories we may replace 
/j«P ^ctp in //ct^ which means this vector is purely spatial with unit norm with respect to g^p. 

Let us look at short wavelength (|) perturbations of RAQUAL in a local Lorentz frame of the 
physical metric gap; the frame is oriented so that 

H'' = e^{0,\,0,0}, (9.2) 

where the e'^ is needed because g^f^H'^H^ = e^'^g^pH'^H^ = e^'^. For g^P we must evidently write in 
Eq. (9J_) e^'''ri"P. Making a WKB approximation, which amounts to replacing 3oi9p8(|) — ^a^p5(|) 



{ka is the wavevector), and discarding the much smaller • • • term, we get the dispersion relation 

(0 = -A:, = [(l+2^)^/ + ^/ + A:,2]V2, (93) 

The group speed Vg = |9(o/9k|'/^ turns out to be anisotropic: 

■(l+2^)V + ^/ + ^; 



(1 +2^)^,2 + ^_^2 +^^2 



1/2 

(9.4) 



In the deep MOND regime [f{y) = ^y^^^], 2^ = 1 while in the high acceleration limit [f{y) ~ 
y], ^ 0. Thus whatever the choice of /, < ^ < 1 over some range of y (acceleration). There 
Vg > 1 unless k^f = 0. This is no artifact of the units chosen; the conformal relation between g„p 
and gap preserves ratios of space and time coordinates, and hence leaves velocities invariant. Now 
light waves travel on light cones of g^p while metric waves do so on null cones of g^p. Since two 
metrics are conformally related, so their null cones coincide, light and metric waves travel with 
unit speed. Thus most (|) waves are superluminal, in violation of the causality principle. It may be 
seen that the blame rests on the vector which, of course, cannot be avoided within this structure 
of RAQUAL. 

Now turning to TeYe S we rewrite Eq. (p?T|) with help of Eq. (|5.3[) as 



[e-2^|«P _ (2 - + 2^«//P] dadp + • • • = 0. (9.5) 

In the local Lorentz frame of g^p whose time and x axes are oriented with and respectively, 
we must put g"P rj^P. Normalizing It" and //" with respect to g^p means taking 

//« = e-'^{0, 1,0,0}; i^« = ^'^{0, 1,0,0}. (9.6) 
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Again making a WKB approximation gives the new dispersion relation 
From this we get the group velocity 

:1+2^)2k,2+k/+K,2- 



-2ij) 



V2 



(9.8) 



The ^ parameter is computed from Eq. ( |5.13[ ): 

^(a/) = (a/-1)(a'-2)/(3a/2-6/j + 4) (9.9) 

from which follows that in the range < /u < I corresponding to y > (recall, in the static case 
y = ^£^|V(j)p), < ^ < 2 with the lower bound corresponding to the Newtonian regime while the 
upper bound is approached in the extreme MOND regime. Consequently, in the deep MOND 
regime, Vg < e^^'^ with equality for ky = k^ = 0. In the Newtonian regime Vg = for all k. 

Finally, in the intermediate regime 2~^l'^e^'^^ < Vg < (1 +2^)'/^2^'/^e"^'l', with lower and upper 
equality for = and ky = k^ = 0, respectively. Since cosmological models with (|) positive and 
small throughout can be constructed in TeYe S, and static systems embedded in them have a (|) with 
the same characteristics (|) waves can propagate causally. 

Causal propagation is maintained in a cosmological background. In addition, it is maintained 
for metric and vector perturbations both on static backgrounds and in cosmological ones [Q]. Thus 
in contrast with RAQUAL, TeYe S is consistent with causality. 



10. Conclusions 

Many a lecture on galaxy phenomenology pointed out that there does not exist a relativistic 
formulation of MOND, and that consequently MOND does not have to be taken seriously as new 
physics. More than anything, the present talk makes it clear that there is nothing that prevents a 
covariant formulation of MOND at the phenomenological level. TeYe S not only reproduces the 
MOND paradigm in the appropriate regime. It also passes the elementary post-Newtonian solar 
system tests, and removes a glaring observational problem with gravitational lensing by extragalac- 
tic systems. And it resolves the problems of superluminal propagation that accrued to aspirants to 
the title of relativistic MOND. But some problems, such as the failure to achieve a perfect Newto- 
nian limit in the outer solar system exist. There remains a large labor to assess how these may be 
fixed, and to extract consequences of TeYe S for the study of cosmological perturbations, gravita- 
tional wave astronomy, binary pulsar timing, and post-Newtonian tests regarding preferred frame 
effects, to name the most obvious. 
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